Exact results on the dynamics of multi-component Bose-Einstein condensate 
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We study the time-evolution of the two dimensional multi- 
component Bose-Einstein condensate in an external harmonic 
trap with arbitrary time-dependent frequency. We show ana- 
lytically that the time-evolution of the total mean-square ra- 
dius of the wave-packet is determined in terms of the same 
solvable equation as in the case of a single-component con- 
densate. The dynamics of the total mean-square radius is 
also the same for the rotating as well as the non-rotating 
multi-component condensate. We determine the criteria for 
the collapse of the condensate at a finite time. Generalizing 
our previous work on a single-component condensate, we show 
explosion-implosion duality in the multi-component conden- 
sate. 

PACS numbers: 03.75.Fi, 05.45.Yv, 11.15.-q, 03.65.Ge 



The successful creation and observation of Bose- 
Einstein condensation(BEC) in dilute alkali atoms have 
opened up a plethora of new possibilities to test, oth- 
erwise intractable, many-body quantum phenomenon in 
the laboratory [Q. The Gross-Pitaevskii equation(GPE), 
the mean- field description of the BEG, is successful 
enough in explaining most of the observed results as well 
as predicting new phenomenon. The methods involved 
in studying the GPE are mainly numerical and/or ap- 
proximate: perturbative and variational. The exact and 
analytical results of a nonlinear equation, if known, not 
only act as a guide to determine the validity of different 
approximate and numerical methods; they also give rise 
to new, counter-intuitive results in some cases. Unfortu- 
nately, no exact solution of GPE is known except for in 
one dimension. 

The two dimensional GPE, like its counterparts in 
higher dimensions, is not exactly solvable. However, due 
to an underlying dynamical 0(2, 1) symmetry ||^, the 
time-evolution of certain moments related to the two di- 
mensional GPE can be described exactly 1^. This result 
is valid even if the condensate is considered in a time- 
dependent harmonic trap. This leads to the prediction of 
explosion-implosion duality Q and extended parametric 
resonance [^,|[ in the two dimensional BEG. Both of these 
phenomenon are universal for any non-relativistic the- 
ory having dynamical 0{2, 1) symmetry [^,^,^. Interest- 
ingly enough, apart from the two dimensional BEG, the 
same explosion-implosion duality can also be observed 
in supernova explosion and in laser induced implosion in 
plasma [0,1). This shows the importance of exact meth- 
ods, based on an underlying symmetry, in relating diverse 
areas of physics such as the BEG and the supernova ex- 



plosion. 

The results described above are for a single-component 
condensate, where the spin degree's of freedom have been 
frozen though the use of a magnetic trap. Recently, the 
spinor condensate with independent spin degree's of free- 
dom has also been created and observed in the labora- 
tory Similarly,the two-component condensate, where 
two different hyperfine states of the same atomic species 
are condensed simultaneously, has also been experimen- 
tally realized [|lO|. The spinor condensate |ll-lj] and the 
two-component condensate [ p^ have a very rich struc- 
ture compared to the single-component condensate. This 
is manifested in the existence of topological defects like 
skyrmion, domain- wall, vortices and Alice string in these 
condensates [ p^ . 

The purpose of this note is to extend the studies 
of Refs. on the two dimensional single-component 
BEG to the two dimensional multi-component BEG. 
The experimentally realizable two-component and the 
spinor condensate can be obtained as special cases of 
this general multi-component BEG. We study the exact 
time-evolution of the second moment of the two dimen- 
sional multi-component condensate in an arbitrary time- 
dependent harmonic trap. This particular second mo- 
ment can be identified as the total mean-square radius of 
the condensate. We show that the dynamics of the second 
moment is determined by the same solvable equation, as 
in the case of a single-component condensate. No matter 
how many components are there, or how they interact 
among themselves, or even whether they are rotating or 
non-rotating, the dynamics of the total mean-square ra- 
dius is universally determined by the same equation. The 
detail information on the system is encoded, through the 
Hamiltonian, into a constant of motion appearing in this 
universal equation. Thus, the dynamics of the system 
can be studied in terms of the same set of initial condi- 
tions for any number of components. We determine the 
criteria for the collapse of the condensate of this system. 
We also show that the multi-component BEG, in its full 
generality, exhibits an explosion-implosion duality and 
extended parametric resonance for special choices of the 
time-dependence of the trap. All these results are exact 
and analytical. 

Gonsider the following Lagrangian in 2 + 1 dimensions. 
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where n is the total number of components. The cou- 
pUng constants gabcd are related to the s-wave scatter- 
ing length matrix. The possible values of gabcd, and 
hence of the scattering length matrix, may be constrained 
by symmetry requirements. For example, the special 
case of a two-component condensate can be obtained 
by choosing n = 2 and Qabcd = ^{Sac^bd + SadSbc)gab 
so that the system has a global U{1)'^ symmetry. A 
phase-separation occurs for such a system if all the 
scattering lengths are positive and satisfy the inequal- 
ity, 512 = 921 > 511522 Similarly, the spin one 
spinor condensate can be obtained by choosing n — 3 

and gabcd = \ [giKc^bd + 52 Y.Jy^a)ac[Sa\i + [o. 6)], 

where S^'s are three spin-matrices. A positive g2 defines 
an anti-ferromagnetic regime, while the ferromagnetic 
regime is characterized by a negative g2 ■ It is known that 
the ferromagnetic or the anti-ferromagnetic nature of the 
interaction plays an important role to characterize differ- 
ent properties of the condensate [jlll,^. Both the phe- 
nomenon of the phase separation and the ferromagnetic 
or anti-fcrromagnctic nature of the ground state are spe- 
cific to multi-component condensate for n > 2. Further, 
note that we have additional terms describing the inter- 
action among different components as we go from the 
single-component to the two-component, to spinor and 
to the general multi-component condensate described by 
(|l|). However, to our surprise, the dynamics of the to- 
tal mean-square radius is independent of such variation 
in the inter-component interaction and universally deter- 
mined by the same solvable equation as in the case of a 
single-component (n = 1) condensate. Consequently, the 
criteria for the collapse of the condensate at a finite time 
is also the same for any n-component condensate. For the 
very special case of an additional global ?7(1)" symmetry 
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in (|l|) , such a result has been obtained previously in Ref. 

We remark that our results are much more general. 
Moreover, the known results are reproduced in a very el- 
egant way. We will consider only the most general form 
of £ from now onwards, since our result is independent 
of particular details of the interaction. 

All the coupling constants gabcd have the inverse-mass 
dimension in the natural units with c — h — I. This 
allows to have a scale and conformally invariant theory. 
The action A = J drd^ rC is invariant under the following 
time-dependent transformations |l8| - [24| , 

dT{t) 
dt ' 



r^rh=T(i) 2r, t i = <(t), f{t) = 



i/'a(T,r) ipa{t,rh) = T-iexp \^-im—rlj -0a(T,r), (2) 
with the scale-factor r given by, 

rW^^rl' «'^-/37=l- (3) 
71 + 

Note that all the components of the order parameter are 
multiplied by the same time-dependent scale-factor and 



the phase in the symmetry-transformation above. One 
might naively think that the requirement of the iden- 
tical phase-factors for all the components of the order 
parameter is due to the interaction term. However, this 
is not the case. Even if we consider the free theory, i.e. 
gabcd ~ for all values of the indices, the requirement of 
the identical phases in is essential in order it to be a 
symmetry transformation. This is precisely because the 
transformation of the scalar-fields is coupled with that 
of the space-time coordinates. If we choose some special 
values for the coupling constants gabcd such that the La- 
grangian has an internal global symmetry, say for exam- 
ple SU{n), we certainly have the freedom of varying the 
phase-factors up to a global SU{n) rotation. However, 
such an additional internal symmetry are completely de- 
coupled from the symmetry transformations described in 
Eqs. (^) and (||), and do not have any effect on our re- 
sults. 

Let us now introduce two moments Ii and I2 in terms 
of the density p and the current j as, 

a 



jd?v p, 
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We are dealing with a conservative system and the total 
number of particles N{t) = J d^rp is a constant of mo- 
tion. The global U{1) symmetry of C can be enlarged 
to for certain special choices of gabcd- The total 

number of particles for each species are conserved sepa- 
rately for this case. However, as emphasized earlier, such 
an additional internal symmetry do not have any signif- 
icant effect on our results. Thus, only the conservation 
of the total number of particles N is important for our 
study. The moment Ii is the sum of the mean-square 
radii corresponding to each and every components. This 
moment can be interpreted as the square of the width 
of the wave-packet for the single-component condensate, 
when confined in an external harmonic trap How- 
ever, for the multi-component case, the moment Ii can 
not be identified as the total width of the wave-packet. 
As emphasized in our previous work the moment Ii 
has been used extensively in the analysis of the non- linear 
Schrddinger equation( NLSE) |,|j2|-|2|, EEC [|§ and 
in optics . The dynamics of Ii , when the system (|l|) is 
immersed in an external time-dependent harmonic trap, 
is the central subject of the investigation of this letter. 
We show that the dynamics of Ii is universally deter- 
mined by the same solvable equation, as in the case of a 
single-component BEC. 
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Particular choices of r(t) = t + f3, a^t, and j:^, corre- 
spond to time translation, dilation and special conformal 
transformation. The corresponding generators of these 
transformations, the Hamiltonian H, the dilatation gen- 
erator D and the conformal generator K are. 



L a abed 



D = tH -U 



(5) 



K = -T^H + 2tD + h. 

These generators close under the algebra, 

[H, D] = iH, [H, K] = 2iD, [K, D] = -iK, (6) 

if we promote the fields ipa 's to the operators ipa with the 
following bosonic commutation relations among them- 
selves, 



Mr),r,{r')\ =5ab <^(r-r'), 
^„(r),^b(r')l = k(r),^b(r') 



0. 



(7) 



The algebra given by Eq. (^) defines a conformal group, 
which is isomorphic to the group 0(2, 1) Thus, the 
system has a dynamical 0(2, 1) symmetry with the 
interpretation of the fields ipa & as the operators -00 sat- 
isfying (^. In this note, we will be considering only the 
fields V'a's, not the operators ^a- We do not make use of 
the relations (|^) or the algebra given by Eq. (^) in our 
subsequent discussions; what is required for our study is 
the conserved Noether charges if, D and K. We just 
mention, in passing, that the results described in this 
note are valid for any non-relativistic theory with a dy- 
namical 0{2, 1) symmetry. 

The generators H, D and K are constant in time and 
lead to the following equations. 



dH dli dl2 

— = 0, -7^ = 2/2, -j^ = H. 
dT dr ar 



(8) 



For time independent solutions, both Ii and I2 do not 
depend on r. As a consequence, the energy of the static 
solutions of H vanishes. This is also the case for the 
single-component BEC in 2 + 1 dimensions. Even though 
there are extra terms due to inter-component interaction 
in the case of multi-component BEC, the vanishing of 
the energy is a universal consequence of the underlying 
0(2, 1) symmetry. The second equation of (||) shows that 
the moment I2 is proportional to the time- variation of the 
moment Ii. Recalling that the moment Ii is identified 
as the total mean-square radius of the condensate, the 
moment I2 can be related to the speed of the growth of 
the condensate. This interpretation is also evident in the 
definition of I2 in Eq. (^) after decomposing the current 
j as a product of the density p and the velocity. 



Defining X = v7i , it is easy to find a decoupled equa- 
tion for X from (pf), 



cPX_ _ Q_ 

dT^ " X3 ' 



= hH-ll ^=0. (9) 



The constant of motion Q is the Casimir operator of 
the 0(2, 1) symmetry. Note that the information on the 
Hamiltonian H is solely contained in Q. Thus, the effect 
of the interaction, say for example the strongly repul- 
sive or attractive inter-component and intra-component 
interaction, will be manifested through initial conditions 
on H. Eq. (^ can be interpreted as the equation of mo- 
tion of a particle moving in an inverse-square potential. 
Interestingly enough, this system also has a dynamical 
0(2, 1) symmetry. This reduced system of a particle in 
an inverse-square potential is a well-studied problem and 
the solution is given by jTsI, 



(10) 



where a and b are integration constants. Although any 
exact solution of the equation of motion of the action A 
is not known, it is surprising to note how the exact time- 
dependence of the moment Ii can be obtained easily us- 
ing the underlying symmetry. We would like to stress 
that we are able to determine the dynamics of the total 
mean-square radius of the condensate only. The dynam- 
ics of the individual mean-square radii associated with 
each components can not be obtained using our method 
even when there is an additional U{1)^ symmetry in the 
system or there is no inter-component interaction. 

The criteria for the collapse of the condensate at a fi- 
nite and real time r* is Q < 0. In particular, the moment 
X^ vanishes at a finite time r*, 



-a6± 



(11) 



which is real if Q < 0. Note that we have the freedom 
of making t* either positive or negative by choosing ap- 
propriate values for the integration constants a and b. 
With the interpretation of (^ as a particle moving in an 
inverse-square potential, the collapse of the condensate 
can be understood as the fall of the particle to the center 
for attractive interaction. Recall that the moment Ii is 
semi-positive definite by definition. Thus, the exact ex- 
pression for Q implies that the condensate collapses for 
any initial condition if < 0. On the other hand, if 
H > 0, the condition for the collapse is given by. 



dh 
dr 



T = 



< -2v//i |r=0 I H 



(12) 



We have used the second equation of (^) in the ex- 
act expression for Q in deriving the above equation. 
As far as we are aware of, this is the first instance in 
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the literature where the criteria for the collapse of the 
condensate of the most general two-dimensional multi- 
component NLSE with cubic nonlinearity is given. The 
criteria is independent of the total number of components 
n and any additional global internal symmetry. Thus, the 
well-known results on the single-component [ p5[ | and the 
multi-component NLSE in two dimensions are easily 
reproduced from our general result. 

Consider the following time-dependent transformation. 



T-^t^t{T), T{t) 

r ^ rh = T{ty^ 



~~ dt ' 
cosf{t) sinfit) 
—sinf{t) cosf{t) 



ipa{T,r) ipait,rh) = iUxp y-'^'^^'^hj ^a(T,r), (13) 

with arbitrary r(i) and f{t). Note that this transfor- 
mation can be obtained by first using the transformation 
(^ and then a time-dependent rotation around the z-axis 
with a time-dependent angle f{t). For arbitrary r(i) and 
f{t), the transformation ( [T^ ) is not a symmetry trans- 
formation of the action A, instead, it maps A to a new 
action Ah = j dtd^vLh- The new Lagrangian Lh reads 
as, 

abed 

- J2 (lmu;it)rl\ |' + , (14) 



where the z-component of the angular momentum, Lz — 
—ivh X Vft and the time-dependent frequency Lij{t) of the 
harmonic trap is determined as. 



b + uj{t)b = {), b{t) = T- 



(15) 



The Lagrangian Ch is that of a rotating multi-component 
BEC in an arbitrary time-dependent harmonic well. Note 
that the external harmonic potentials are identical for 
all the components of the condensate. This is not by a 
choice. In fact, we do not have the freedom of generating 
different harmonic potentials for different components us- 
ing the transformation in Eq. (p^). This is even true for 
the free theory, i.e. gabcd = 0. The reason is that the 
transformation of the scalar fields is coupled with that 
of the space-time coordinates. Consequently, unphysical 
and unwanted terms will be generated in the new La- 
grangian Ch unless all the components of the condensate 
transform identically. 

The solutions of A and Ah are related to each other 
through the transformations in Eq. (^) with rit) deter- 
mined for a specific trap-frequency by the Eq. (p^). The 
scale-factor r(t) can obviously be exactly determined for 



a large class ofuj{t). However, the exact solutions are not 
known for either A or Ah- This is a major problem in 
making use of the mapping relating A to Ah and the vice 
versa. However, note that the dynamics of the moment 
Ii is uniquely determined by Eq. ( pO| ) independent of 
whether any exact solution of A is known or not. Thus, 
the transformation ( p^ can be used to find the dynam- 
ics of the moment h^h — J2a I '^'^^h'rf^l ipa f from Ji. In 
particular, they are related to each other by the relation. 



b{t)X{T{t)), 



(16) 



where b{t) and T{t) are determined from the Eq. (|l|). 
Thus, even though the exact solution of the equation of 
motion of Ah is not known, the dynamics of Xh can be 
described exactly. 

An alternative, but, useful expression for the Xh can 
be determined from the following equation M, 



fXh 
dt^ 



+ Lo{t)Xh = Qh = IiMHh - Ilh, (17) 



where Qh is a constant of motion. Both Hh and l2,h have 
the same expressions as in H and I2, respectively, with 
(r, r,V'a) replaced by {t,rh,ipa)- Note that the Eq. ( |T7| ) 
can also be interpreted as describing the motion of a clas- 
sical particle in a combined harmonic and inverse-square 
potential. The particle falls to the center for an attrac- 
tive {Qh < ) inverse-square potential, independent of 
the time-dependence of the harmonic trap. This implies 
that the condensate collapses at a finite time for Qh < 0. 
Analyzing the exact expression for Qh further, we find 
that the condensate collapses for any initial condition if 
Hh < 0. For Hh > 0, the condition for the collapse is 
given by. 



dlij 
dt 



< 



Hh 



\t=o, 



(18) 



where the relation ii_h = 2/2, /i, valid for the system 
described by the action Ah has been used. As far as 
we are aware of, this is for the first time in the litera- 
ture that a criteria for the collapse of the condensate in 
the most general two dimensional multi-component GPE 
with cubic nonlinearity and an arbitrary time-dependent 
harmonic trap is given. Note that the criteria is inde- 
pendent of the total number of components n and any 
additional internal global symmetry. The known results 
for the single-component and the multi-component 
jlTf GPE with time-independent harmonic trap in 2 + 1 
dimensions are easily reproduced from this very general 
result. Further, the criteria for the collapse in the sys- 
tem without or with the harmonic trap is also identical, 
except that there is no equality sign in (|l8|) for the later 
case [ compare with Eq. (p^)] . This is precisely because, 
for Q, Qh = 0, Eqs. (^) and ( p7[ ) describe the dynamics of 
a free particle and that of a particle in a time-dependent 
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harmonic trap, respectively. Thus, nothing can be said 
conclusively on the dynamical (in)stability for the later 
case, unless the time-dependence of the frequency of the 
trap is explicitly specified. For a time-independent trap, 
the equality sign is recovered in Eq. (jl^); and of course, 
the known result |l^,|2^ is identically reproduced. 

We have shown that the criteria for the collapse of the 
condensate in a 2 -I- 1 dimensional system governed by 
the Lagrangian Ch is independent of the total number of 
components n. It is known | |2^ that the same criteria 
for the collapse of the condensate is also valid for the 
Lagrangian Ch in dimensions d > 2 -I- 1 with n = 1, 
Lv{t) = Wo = a constant and / = 0. So, the criteria 
for the collapse is independent of the underlying 0(2, 1) 
symmetry, which the cubic NLSE has only in d = 2 + 
1. The dynamical 0(2, 1) symmetry only helps us in 
deriving the exact result in a much more simpler and 
elegant way. Based on this observation, we conjecture 
that the criteria for the collapse of the condensate of Ch 
in dimensions d > 2-1-1, with aj(t) = loq and / = 0, 
is independent of the total number of components n and 
the criteria is the same as stated in this note for c? = 2-1-1. 
Note that the physically interesting case of d = 3 -I- 1 is 
also included in our conjecture. 

The solution for Xh is given by, 



square radius. In particular, the dynamics of the moment 
Ii is determined from the same solvable equation as in 
the case of a single-component BEC with all the infor- 
mation about the Hamiltonian encoded into the constant 
of motion Qh- Apart from its relevance to the ongoing 
experiments on BEC, the importance of these results lie 
in its universality. No matter how many components are 
there, or how they interact among themselves, or even 
whether they are rotating or non-rotating, the dynamics 
of the total mean-square radius is universally determined 
by the same equation. This is indeed a counter-intuitive 
result and may be realized in the laboratory in near fu- 
ture. 
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where u{t) and v{t) arc two independent solutions of Eq. 
(|l|) satisfying u{ta) = Xh{to), u{to) = Xh{to), v{to) = 0, 
and w(to) 7^ 0. The above solution is valid for arbitrary 
Qh '■ positive, negative or zero. We will be considering 
the case Q/i > from now onwards, since we have al- 
ready argued that the condensate collapses for Qh < 0. 
We have obtained the same expressions ([l^ ) and (|9|) in 
Refs. for the dynamics of the width of the wave- 
packet of a single-component condensate in 2 -I- 1 dimen- 
sions. So, the results of the Refs. ||,^ are equally valid 
for the general multi-component condensate in 2 + 1 di- 
mensions with the moment Ii — Xf^ identified as the total 
mean-square radius. In particular, 

(a) the system described by Ch has an explosion- 
implosion duality for fit) = and either Lo{t) — or 
io{t)^t-\ 

(b) the condensate exhibits extended parametric reso- 
nance for a periodic uj[t) and arbitrary /(t), 

(c) the dynamic (in) stability of the system is independent 
oi f(t), i.e., same for both rotating as well as non-rotating 
BEC. 

We refer the readers to Refs. for further details. 

Finally, we conclude with the following comment. The 
results presented in this note for the multi-component 
BEC are a generalization of what already has been known 
for the single-component BEC in two dimensions. The re- 
sults obtained in both these cases are also identical with 
the identification of the moment Ii as the total mean 



[1] 
[2] 



F. Dalfovo, S. Giorgini, L. P. Pita evskii and S. String ari, 
Rev. Mod. Phys. 71 463 (1999) (Jcond- mat /980603f ) . 
We refer to the papers in Refs. [[isf and [|19|| for the appli- 
cation of 0(2, 1) symmetry in physics. The models con- 
sidered in these papers also have a dynamical 0{2, 1) 
symmetry. 



[3] P. K. Ghosh Phys. Rev. A65, 012103 (2002) ( |cond 



w 

[5] 



mat/0102488 



:ond-mat/0109073| . 



P. K. Ghosh, 

J. J. Garcia- Ripoll an d V. M. Pcrcz-Ga rcia, Phys. Rev. 



Lett. 83, 1715 (1999); patt-sol/9904006 



mat/0012188 



[61 T. K. Ghosh , Phys. Lett. A 285, 222 (2001 ) ( |cond 

[8 



I. O'C Drury and J. T. Mendonc a, Physics of Plasmas 7, 
5148 (2000) ( |astro-ph/0003385| ). 

K. Walter, Supernova Hydrodynamics Up Close: 
Science and Technology Review, Jan'/Feb' 2000 ( 
www.llnl.gov/str/BRemington.html); L Hachisu et. al., 
Astrophysical Journal 368, L27 (1991); H. Sakagami and 
K. Nishihara, Physics of Fluids B 2, 2715 (1990). 
[9] J. Stenger, S. Inouye, D. M. Stamper-Kurn, H.-J. Mies- 
ner, A. P. Chikkatur and W. Ketterle, Nature 369, 345 
(1998); H.-J. Miesner, D. M. Stamper-Kurn, J. Stenger, 
S. Inouye, A. P. Chikkatur and W. Ketterle, Phys. Rev. 
Lett. 82, 2228 (1999); D. M. Stamper-Kurn, H.-J. Mies- 
ner, A. P. Chikkatur, S. Inouye, J. Stenger and W. Ket- 
terle, Phys. Rev. Lett. 83, 661 (1999). 
[10] C. J. Myatt, E. A. Burt, R. W. Christ, E. A. CorneU 
and C. E. Wieman, Phys. Rev. Lett. 78, 586 (1997); M. 
R. Matthews, D. S. Hall, D. S. Jin, J. R. Ensher, C. E. 
Wieman, E. A. Cornell, F. Dalfovo, C. Minniti and S. 



5 



[11] 

[12] 
[13] 

[14] 
[15] 



[16] 



Stringari, Phys. Rev. Lett. 81, 243 (1998); D. S. Hall, 

M. R. Matthews, J. R. Ensher, C. E. Wieman and E. A. 

Cornell, Phys. Rev. Lett. 81, 1539 (1998). 

T. Ohmi and K. Machida, J. Phys. Soc. Jpn. 67, 1822 

(1998). 

T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998). 
D. M. Stamper-Kurn and W. Ketterle, 



cond-mat/0005001 



[17] 
[18] 
[19] 
[20] 



A. J. Leggett, Rev. Mod. Phys. 73, 307 (2001). 
T.-L. Ho and V. B. Shenoy, Phys. Rev. Lett. 77, 3276 
(1996); B. D. Esry, C. H. Greene, J. P. Burke, Jr. and J. 
L. Bohn, Phys. Rev. Lett. 78, 3594 (1997). 
S.-K. Yip, Phys. Rev. Lett. 83, 4677 (1999); U. Leon- 
hardt and G. E. Volovik, JETP Lett. 72, 46 (2000); J. 
Ruostekoski and J. R. Anglin, Phys. Rev. Lett. 86, 3934 
(2001); H. T. C. Stoof, E. Vliegen and U. Al Khawaja, 
Phys. Rev. Lett. 87, 120407 (2001); U. Al Khawaja and 
H. Stoof, Nature 411, 918 (2001); K.-P. Marzhn, W. 
Zhang and B. C. Sanders, Phys. Rev. A 62, 013602 
(2000). 

T. Tsurumi and M. Wadati, J. Phys. Soc. Japan 67, 93 
(1998). 

V. de Alfaro, S. Fubini and G. Furlan, Nuvo Cimento 
A34, 569 (1976). 

R. Jackiw, Ann. Phys. (N.Y.) 129, 183(1980); 201, 83 
(1990). 

C. R. Hagen, Phys. Rev. D5, 377 (1972). 



U. Niederer, Helv. Phys. Acta 45, 802 (1972); 46, 191 
(1973); 47, 119, 167 (1974); 51, 220 (1978). 
E. A. Kuznetsov and S. K. Turitsyn, Phys. Lett. A112, 
273 (1985). 

S. Takagi, Prog. Theor. Phys. 84, 1019 (1990); 85, 463 
(1991); 85, 723 (1991). 



[21] 
[22] 
[23] 

[24] Pijush K. G hosh, Phys. Rev. D 53, 2248 (1996) ( |hep 

[25; 



th/9503199 



V. E. Zakharov, Sov. Phys. JETP 35, 908 (1972); V. E. 
Zakharov and V. S. Synakh, Sov. Phys. JETP 41, 465 
(1975); M. L Weinstein, Commun. Math. Phys. 87, 567 
(1983). 

[26] T. Tsurumi and M. Wadati, J. Phys. Soc. Japan 66, 3031 
(1997); ibid 66, 3035 (1997); M. Wadati and T. Tsurumi, 
Phys. Lett. A 247, 287 (1998); T. Tsurumi, H. Morise 
a nd M. Wadati, Int Jour. Mod. Phys. B 14, 655 (2000) 
( |cond-mat/991247C| ). 

[27] A. V. Rybin, G. G. Varzugin, M. Lindberg, J. Timonen 
and R. K. Bulloug h, Phys. Rev. E 62, 6224 (2000) ( 
cond-mat/0001059 ). 

[28] S. Stringari, Phys. Rev. Lett. 77, 2360 (1996); V. M. 
Perez-Garcia et. al., Phys. Rev. Lett. 77, 5320 (1996). 

[29] S. N. Vlasov, V. A. Petrischev and V. L Talanov, Radio- 
phys. Quantum Electron. 14, 1062 (1971); V. M. Perez- 
Garcia et. al., J. Opt. B : Quantum Semiclass. Opt. 2, 
353358 (2000). 



6 



